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Quantum transport in a lattice is distinct from its counterpart in continuum media. Even a free 
wave packet travels differently in a lattice than in the continuum. We describe quantum scattering in 
a one dimensional lattice using three different formulations and illustrate characteristics of quantum 
transport such as resonant transmission. We demonstrate the real time propagation of a wave packet 
and its phase shift due to impurity configurations. Spin-flip scattering is also taken into account 
in a spin chain system. We show how individual spins in the chain evolve as a result of a spin-flip 
interaction between an incoming electron and a spin chain. 

PACS numbers: 72.10.-d, 73.50.Bk, 73.21.Hb 
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I. INTRODUCTION 

In the last few decades the advent of pump-probe op- 
tical methodsiiSi^ and similar measurement techniques 
has stimulated interest in time-dependent phenomena 
in physical systems^ For example, questions concern- 
ing the details of magnetization reversal in ferromag- 
netic thin films can be addressed experimentally^, and 
theoreticallji^. Another example is polaron formation in 
a semiconductor, where both experimental and theoret- 
ical work are just starting.^ While much of the work on 
magnetization reversal has utilized a classical or at most 
semi-classical description(2i2ii24iiii24i^ more recent work 
has focused on a fully quantum mechanical description J4. 
Such a microscopic description is expected to be neces- 
sary and insightful for small (i.e. quantum dot) systems. 

Insofar as many of these phenomena occur in the solid 
state, the underlying lattice structure may play an im- 
portant role. Concerning the theoretical description of 
magnetization reversal, calculations to date have either 
focused on simple models^SiSiiSiiiii^ or have tried to uti- 
lize realistic transport equations with band structures 
relevant to the materials of interest ii& In this work we 
wish to solve simple scattering models based on tight- 
binding band structures. While we will consider mainly 
scalar potential scattering, the formalism is extendable 
to spin-flip scattering, which will be discussed. In par- 
ticular, an incoming electron can scatter off a ferromag- 
netic thin film modelled by a Heisenberg Hamiltonian, 
and we can monitor the real time reaction of the magne- 
tization to the onslaught of electrons with a completely 
quantum mechanical description. We also want to utilize 
a framework that is amenable to numerical calculation. 
By this we mean the following: as interactions are in- 
troduced, problems will become formidable by analyti- 
cal means and large scale computation will be required. 
Most often this means Monte Carlo methods (at least 
this is so in equilibrium and linear response calculations 
so far) which are often well-suited to simple lattices. For 
these reasons we believe it is beneficial to have a lattice- 
based framework to tackle non-equilibrium phenomena 



in solid state systems. 

We begin with a description of the non-interacting elec- 
tron, where already a novel property emerges due to the 
lattice: the degree of spreading of a propagating wave 
packet can be controlled by judicious choice of the elec- 
tron energy. This is always true in one dimension and 
has more limited validity in higher dimensions. Next we 
solve the scattering problem for a rectangular potential 
barrier and other, simple barriers (or wells) that enter 
into impurity problems. We outline the methodology to 
solve the problem numerically, and use some illustrative 
examples to demonstrate the accuracy and efficiency of 
these calculations. Finally, we discuss spin ffip scattering 
in a spin chain. 

This paper is organized as follows. In Section. II, we 
describe wave packet transport in a lattice without scat- 
tering, and a possible solitonic behavior which is impossi- 
ble in the continuum limit. In order to illustrate quantum 
scattering off impurities in a lattice, we use a quantum 
mechanical approach and the transfer matrix formalism, 
and compare the two methods using an example in Sec- 
tion III. A direct diagonalization of the Hamiltonian will 
illustrate the time evolution of a wave packet in a lattice 
with embedded impurities in Section IV. In Section V, 
we briefly explain a procedure to study spin flip scatter- 
ing on a lattice. In section VI we summarize our results 
and outline possible future directions. 



II. WAVE PACKET TRANSPORT 

In order to understand the differences between calcu- 
lations of quantum mechanical phenomena in a lattice 
from their counterpart in the continuum limit, it is ap- 
propriate to begin with transport of a free wave packet. 
As is normally done in textbooks^ one can introduce a 
free Gaussian wave packet in the continuum limit: 
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where Xq and kg are the mean position and momentum, 
respectively, of the wave packet, and a is the position 
uncertainty associated with the wave packet. To see the 
time evohition of the wave packet we expand ^^(x, 0) 
in terms of the momentum eigenstates |fc). They are 
the known solutions to the Schrodinger equation in free 
space; hence the time dependent Schrodinger equation is 
readily solved in this basis. The result isiii 
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(2) 
where Eq — k'^/2vi and wq — ko/m are the average en- 
ergy and particle velocity, respectively. Since (x) — vat 
and (x^) = (vot)^ + [a^ + {t/2ma)^'\ , the position un- 
certainty Ax defined as (Ax)^ — (x^) — (x)^ becomes 
Ax = y^a"^ + 1'^ /{2ma)^. Similarly, we obtain the mo- 
mentum uncertainty Afc — l/{2a). The uncertainty re- 
lation is therefore 



Ax • AA: 



1 + 



t 



2ma^ 



(3) 



This means that the uncertainty relation increases as a 
function of time. Also note that the relation does not 



depend on the mean momentum of the wave packet; as we 
will see, this is true only for a parabolic energy dispersion. 

In a one dimensional lattice described by a nearest 
neighbor tight-binding model, Ek — —2tocos{ka), where 
to is a hopping amplitude to the nearest neighbor site 
and a is a lattice constant. Hereafter we set a = 1 and 
use it as the unit of length. The position is now discrete 
and represented by Xi, with i the lattice site label. Upon 
expanding in terms of the momentum eigenstates in a 
box with periodic boundary conditions, one obtains 



■^{Xi,t) 
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dk J'^^i^°'^ (^^^of -i^kt /^\ 



Note that the integration is over a Brillouin zone, due to 
the discreteness of the lattice; nonetheless, if a is large 
enough, the integration range can be extended from — oo 
to oo without altering the integral. In the same way, 
using a large a expansion for the exponent in Eq. Q), 
and keeping terms up to 0{l/a^), one can convert the 
integral into a Gaussian integral, as in the continuum 
limit. Performing the integration, we obtain 
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^^{x^,t)= —r 
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where we restored more generality (than in Eq. ^J by 
using El and E'/ to refer to the first and second deriva- 
fives of the dispersion, Ek, with respect to momentum 
k, and evaluated at k^. For a quadratic dispersion one 
readily obtains the result Eq. 12 On the other hand, for 
the nearest neighbor model, E'j, = Vkg = 2iosin(/co) and 
E'j! = v'j, — 2iocos(fco), where Vk^ is the group velocity, 
and v'f, is the group velocity dispersion. The expansion 

is vahd as long as i ^ (a/fco) /to- For /cq = 7r/2, the 
validity length I is I = Vkgi ~ a^. If a is on a nanometer 
scale, I is of order l^m. 

We are now able to calculate the uncertainty relation 
for the lattice case with nearest neighbor hopping only, at 
any time i: (x) — 2tosin(/co)i and (x^) = {2to sm{ko)t)'^ + 
o? + (ig cos{ko)tloi) 
Ax = ^a^ +tlcos'^{ko)t'^/a^. 
momentum is the same as in the continuum limit; namely, 
Afc = l/2a. Consequently, the uncertainty relation for 
this case on a lattice is 



The uncertainty in position is then 
The uncertainty in the 



Ax • Afc = - 
2 



t2c0s2(fco)t2 



1 + 



a^ 



(6) 



As one can see from this expression, in general the un- 
certainty never decreases as a function of time; the de- 
gree of increase depends on the mean momentum. How- 
ever, if fco = 7r/2, the uncertainty remains unchanged. In 
other words, the wave packet possesses a solitonic behav- 
ior without showing the seemingly inevitable quantum 
spreading. We will demonstrate this fact numerically 
later. This possibility is actually well known in optics, 
where one seeks a medium with zero group velocity dis- 
persion to minimize lossii^ Nevertheless, this appears to 
be less appreciated in quantum mechanics. 

This result persists in one dimension for any dis- 
persion. That is, one can show that some wave vec- 
tor always exists for which the group velocity disper- 
sion is zero. In higher dimensions the situation is not 
quite as simple. The result remains for nearest neigh- 
bor hopping only. For example in three dimensions we 
have Ek = — 2to [cos(fc2:a) -t- cos(fcj,a) -t- cos(fcza)], and 
one readily obtains a similar result as in one dimen- 
sion. However, when next nearest neighbor hopping is 
included, a little algebra shows that in general one can- 
not achieve conditions for zero group velocity dispersion. 

To obtain this result numerically, one diag- 



Fig. 1 (Kim et al) 




FIG. 1: Time evolution of the wave packet with fco = 7r/2 
(top) and 7r/4 (bottom). The wave packet with fco = 7r/2 
does not spread while the wave packet ko = 7r/2 broadens as 
it moves. 



onalizes the tight-binding Hamiltonian, Hq = 
— ioX^i [^i^C'i+i + Cj^^Ci] *° obtain eigenvalues and 
the corresponding eigenstates. Then one can construct 
a wave packet, and evolve it in time according to the 
time dependent Schrodinger equation (this procedure is 
described in more detail in Section IV). The result is 
plotted in Fig. 1 for three different times with fco = ""/S 
(fco = 7r/4) in the top (bottom) panel, respectively. A 
wave packet initially centered at a;,; = 100 moves to 
Xi ~ 300 in each panel. The uncertainty parameter a is 
set to be 10. The total time elapsed is different in the 
two panels because of the fco difference (in a consistent 
dimensionless time unit, 100 for the top, and 140 for the 
bottom panel). As one can see, the wave packet with 
fco = 7r/2 does not spread and the peak height remains 
unchanged while the wave packet with fco = 7r/4 spreads 
and the height becomes considerably smaller as it moves. 



III. QUANTUM MECHANICAL APPROACH 

In the previous section we examined the free wave 
packet transfer in a lattice. While construction of a 
packet proceeds as in the continuum case, the tight- 
binding dispersion leads to a variety of possible wave 
packet characteristics. In this section we examine the 
consequences for scattering, in one dimension. The elec- 
tronic transfer can be explored based on the quantum me- 
chanical approach using the matching conditions of wave 
functions, or through the transfer matrix formalism, or 
through direct diagonalization of the Hamiltonian. The 
quantum mechanical approach and the transfer matrix 
formalism are formally equivalent to each other and rel- 
atively analytical compared with direct diagonalization, 
which is inherently numerical and the subject of a later 
section. If there are N impurities embedded in a lattice. 



one needs to deal with a (2A^ x 2A^) matrix to deter- 
mine all the relevant coefficients including R and T in 
the quantum mechanical approach. On the other hand, 
the transfer matrix formalism requires manipulation of 
iV (2 X 2) matrices. Thus, when the number of impuri- 
ties is large, say greater than 5, the quantum mechanical 
approach is not as feasible as the transfer matrix formal- 
ism. 

We used the transfer matrix formalism to study sym- 
metries of electron- impurity scattering previously.^^^ In 
this paper, we will mainly use the quantum mechanical 
approach in a lattice and compare the two methods using 
several examples. Recall the situation in the continuum 
limit: consider a Hamiltonian H = — 0^/(2to) + V{x) 
where, for this illustration we use a simple rectangular 
barrier potential: V{x) = VoQ{d/2 — \x\). In order to 
solve the Schrodinger equation, T-iip{x) = Ek^{x) with 
Ek = fc^/(2?7i), we write down the most general wave 
function in the three regions. 



tpLix) -- 
1p{x) = { 1pin{x) 
'^R.ix) -- 



^ikx 



e'"-" + Re 
-- A e*«^ +B e 

rp ikx 



ikx 



-iqx 



for X < -d/2 
for |a;| < d/2 
for X > d/2 



where q is determined hy q — y^2m{Ek — Vq). Note that 
'^Jl,r(x) describes a free electron satisfying TiotpL.Rix) — 
E^l,r{x) for any value of x, where Hq — —dl/{2m). 
Using the matching conditions of ip{x) and dipix)/dx at 
X = ± d/2, we determine the coefficients R, A, B, and 
T. An examplc"'^^ of this case is the one as d ^ and 
Vb ^ oo while dVo is finite; namely, V{x) — dVoS{x). A 

little algebra yields |r|2 = 1/ l + {dVo/vkf and\R\^ = 

1 — |Tp, where the velocity Vk = dEk/dk = k/m. 

In a lattice, a potential is represented by a series of 
"impurities" whose effect is to alter the on-site energy 
wherever an impurity has substituted for the usual atom. 
This is represented by the Hamiltonian, 



H^-toJ2[CtC' 



^+l + c+,,a]+Y,u.c+a, (7) 

iei 



where to is the hopping amplitude as before, C/' cre- 
ates an electron at a site i, and Ui is a scalar poten- 
tial at site i; the set of "impurities" spans a number 
T — {0, 1, 2, • • • ,/}. Note that Ui can be positive (re- 
pulsive), negative (attractive), or zero. In this way one 
can construct any shape potential one wishes (on a scale 
much greater than the lattice spacing). Since we con- 
sider only scalar potentials, and we are treating the one 
electron problem here, we ignore the spin index of the 
electrons. The wave function defined only on the lattice 
sites can now be written as a piecewise function over 1+3 
regions since the number of impurities is / -|- 1. That is, 
lV')=E,Xa^.)C^|0), where 



ip{x. 



i^bixj) = e 
ikj{xj) ^Aj 
ipF/ixj) = T e 



e^ij-^j + R. e-««j^j for j el 



for j > I 



and |0) represents the vacuum, namely, the state with 
an empty lattice. Since we set the lattice constant to 
be unity, the displacement Xi = a ■ i — i. The coeffi- 
cients R, Ai, Bi, and T are to be determined by match- 
ing conditions at z G X. Note that qi will be obtained 
within the same calculation even though we can already 
guess that E^ = Eg. + Ui. What are the matching con- 
ditions ? Like the continuum limit one first demands 
continuity of the wave function at each site. Therefore, 
(0|Cj+|^) = {Q\Cj-\tp), where j-l- (j— ) means just to the 
right (left) of site j. However, the second condition in 
the continuum limit requires continuity of the derivative 
of the wave functions at each site. One can see this di- 
rectly from the Schrodinger equation through an integra- 
tion of the second order differential equation. However, 
the second quantized form of the Hamiltonian written in 
Eq. ^ contains no derivatives, so clearly this procedure 
is not an option. The correct procedure is as follows. ^"^ 
One first writes down the Schrodinger equation projected 
onto each site, (0|Cji/|V') = {Q\CjE\^) for j = 0, 1, • • • , /. 
Then the two conditions, expressed for each site, can be 
written 



-io [iAi + 1) + V^J - 1)] + U,i;(j) = Ei:{j) . (9) 



^(j + 0^ 



i^ii + 0- 



(8) 



As mentioned before, Eq. (jH)) implies the wave function is 
continuous at each site. Eq. (O is the Schrodinger equa- 
tion at each site; however, close inspection shows that 
on the left-hand side the wave function is required from 
two different "pieces" in the domain (assuming that Uj is 
non-zero. But we would like the Schrodinger equation for 
non-interacting electrons to be satisfied, with the same 
eigenvalue, by the wave function on any given "piece" 
even when extended beyond the domain of validity of 
that wave function. For example, for j — Q we demand 
that Vl satisfy the equation, —to [V'l(+1) +'0l(— 1)] — 
EiPl{0). Note that we have used the wave function ■0^ at 
location -1-1, even though it was originally defined only 
for sites or below. Moreover, we require that this equa- 
tion be satisfied with the same eigenvalue, E. Hence, by 
judicious adding and subtracting of a wave function to 
Eq. (jSJ at each impurity site, we arrive at, for j = and 



-iol/^l(l)+ioV'L(l) 
-to4>liI-l)+to4,R{I-l) 



-UoiJLiO) = forj=0 
Ui^Pnil) = forj=/ 



r 



(10) 

(11) 



Similar equations apply for the impurity sites in between. 
These can now be solved for the unknown coefficients 
along with the "matching equations" [Eq. ©]. 

Let us consider a two impurity case as an example. 
Assume two impurities with Uq and Ui are embedded 
at J = and 1, respectively, in a lattice; one needs to 
introduce a wave function as follows: 



^dxj) 






Re 






lP{Xj) = 



V'o(a;j) = Ao e'?«="^- + Bq e-**"^^" 
%Pi{x/) = Ai e*?i^^- + Bi e-*«i^^" 
t^Rix,) = Te'''-^. 



The coefficients such as R and T, and the momenta 
qq, and qi can be determined by solving the continuity 
equations Vl(-1) = V'i(-l), i/'l(0) = MO), MO) = 
MO), ^i(l) = M^), V'2(l) = ^fl(l), ^i(2) = Vfl(2), 
and the Schrodinger equations 



>i(1) + ^l(1) 

-M0) + M{0) 



Uo^l{0) 

Uli^Ril) 





0, 



(12) 



where we set the nearest neighbor hopping amplitude io 
to be unity for simplicity. It is straightforward to show 
that 



T 



2isin(fc) 



2i sin(fc) - (C/o + C/i) - UoUie^'' 



(13) 



R = 



Uo + Uie^"" + UoUie 



ik 



2i sin(fc) - {Uo + Ui) ~ UqUic^^ ' 



(14) 



Using the Schrodinger equations, one can show that 
Ek — Eq. +Ui {i = 1,2). This result will be used to com- 
pare the quantum mechanical approach and the transfer 
matrix formalism. A case of two impurities with same 
potentials (t/o = ^i) has been well studied in a context 
of the random dimer modeli^SiSiSi 

In Fig. 2, we plot jTp for two impurities side by side 
using Eq. l(T^ . Fig. 2(a) to Fig. 2(c) are contour plots 
of \T{E,Ux)\^ for given values of Uo = -1, 1, and 1.8, 
respectively. Note that Fig. 2(a) and Fig. 2(b) illustrate 
a symmetry, \T{E, Uo, C/i)P = \T{^E, -Uo, -C/i)|2, 
which we showed in a previous paper.^-'^ When we in- 
crease Uo further, |Tp becomes suppressed considerably 
as shown in Fig. 2(c). We also plot |Tp vs. E in Fig. 2(d) 
for various values of Ui while C/q is fixed to be 1.0. In 
fact, Fig. 2(d) can be read off from Fig. 2(b) along a line 
of the corresponding value of Ui . In the case of two impu- 
rities, a transmission resonance (|rp = 1) occurs when 
E = Uo = Ui and \E\ < 2. As one can see in Fig. 2(d), 
a transmission resonance (|rp = 1) occurs at E — I for 
Uo = Ui = 1. 
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FIG. 2: [Tp for two impurities using Eq. ^. Fig. 2(a) 
to Fig. 2(c) are contour plots of T{E,Ui)\ for given values 
of Uo = —1, 1, and 1.8, respectively. Fig. 2(d) is \T\^ as a 
function of E for various values of Ui with Uo = 1.0. 



For three impurities with Uo, Ui, and U2 at j — 0, 1 
and 2, respectively, a similar calculation gives the trans- 
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mission amplitude 



T = 



-2ism{k) 



(C/i - Ek) + e^^ [(t/o + U2){Ui - Ek) - 2] + e^^^ [UoU2{Ui - E^) - {Uo + U2)] ' 

I 



(15) 



It is easy to shovi' that Eq. H15|) reduces to Eq. (|13|) as Uo 
or U2 goes to zero. Giri et alm^ studied this case in detail 
and classified the parameter values into five cases, for 
which perfect transmission can be obtained. The first is 
the dimer case, and the second only applies if the hopping 
parameters in the impurity region are different, so here 
we will consider only the remaining three cases. Let us 



briefly describe the three cases using our notation. In 
case (III), Uo ^ Ui = U2 and E = Uo ± 1- This is the 
straightforward extension of the symmetric dimer case to 
a symmetric trimer, but note that now the potential for 
two solutions exists. In case (IV), if Uo + U2 = Ui then 
resonance is obtained at an energy, E = Ui, regardless 
of the values of Uo and C/2. We will note an example 
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FIG. 3: |rp for three impurities using Eq. ^^. Fig. 3(a) to Fig. 3(c) are contour plots while Fig. 3(d) is for \Tf vs. E. 
Fig. 3(a) is for |T(£,f/o)p with !7i = C/2 = 0.6. Fig. 3(b) is also for \T{E,Ua)\^ but with Vx = 0.6 and V2 = 0.8. Fig. 3(c) 
shows |r(_E, [7i)|^ with Uo — U2 — 0.6. Fig. 3(d) is for \Tf vs. E, which can be read off from Fig. 3(a) along a line of the 
corresponding value of Uo- 



below. Case (V) is the most interesting; in this instance 
Uq = U2 and Ui is arbitrary. Then resonance is obtained 
for energies, 



E = Ui + - 



Ui~Uo± ViUi - Uof - 4{Ui/Uo - 2)J , 

(16) 
provided that (C/i — C/q)^ > 4:{Ui/Uo — 2) and of course 
—2<E<2. We plot |Tp for three impurities sitting 
side by side in Fig. 3. Fig. 3(a) to Fig. 3(c) are contour 
plots while Fig. 3(d) is for jTp vs. E. Fig. 3(a) is for 
\T{E, C/o)P with Ui^U2 = 0.6 and illustrates both case 
(III) and case (IV). For example, if Uq = 0.6, \T\^ = 1 
at E — 1.6 and —0.4. This belongs to case (III). On the 
other hand, since Ui — U2, then when Uq = 0, |Tp = 1 
at _E = 0.6, which belongs to case (IV). In Fig. 3(b), 
which is also a contour plot of \T{E, C/o)P but with Ui = 
0.6 and U2 = 0.8, we still have case (IV) when C/q = 
-0.2 and E = 0.6. On the other hand, when Uq = 0.8, 
|rp = 1 at £; = 1.82 and -0.42, which belongs to case 



(V). Fig. 3(c) is primarily an example of case (V) since 
Uq = U2. The condition {Ui - Uof > 4([/i/[/o - 2) 
requires Ui < 1.27 for Uq = 0.6. Since \E\ < 2, one finds 
that there exist two solutions for jTp = 1 if —0.75 < 
Ui < 1.27; however, for -3 < C/i < -0.75 there is only 
one solution as can be observed in Fig. 3(c). Also note 
that Fig. 3(c) contains as a special case the 'spread-out' 
dimer, where the middle of the trimer has no impurity 
(Ui = 0), so the remaining dimer is now separated by 
two lattice spacings. This is a special example of Giri et 
al. 's case (V) , and is also a special case of the arbitrarily 
spread-out dimer considered by us previously^^ . 

Fig. 3(d) is for |Tp vs. E, which can be read off from 
Fig. 3(a) along a line of the corresponding value of Uq. 
Note the two solutions for perfect transmission for the 
symmetric trimer case (dashed blue curve, Uq — 0.6). 
Also shown is a dimer case (dashed green curve), when 
Uq = 0, for which E = 0.6 yields resonant transmission. 
The standard one dimensional plots are quite clear; on 
the other hand the contour plots provide a feel for how 



strongly the transmission remains as parameters vary 
away from the resonant condition. 



IV. TRANSFER MATRIX FORMALISM 

In the transfer matrix formahsm,^® we write the 
Schrodinger equation Eq. (PJ in the matrix form as fol- 
lows: 



where Mj — 



i'j 

1 



M, 



-1 






(17) 







which is a unimodular 



matrix and associated with an impurity at the site j. 
The wave functions i/'l (for i < 1) and ipji (for i > N) 
are i/^l = e'^'"^^ + R e~'^'"^^ and Vfl = T e^^^^ . Using 
the transfer matrix formalism, one can express the coef- 
ficients R and T in terms of fc, C/i, and E as follows: 



r 

iT 



P 



i(l-7^) 



(18) 



where P = S-'^MS with S 



cos(fc) sin(fc) 



M = MnMn-i---Mi. 
obtainSi 



T 

n 



1 

Solving Eq. (|18|l . one can 



and 



2i 



I [Pii + P22) + P12 - P21 
P12 + P21 - i {Pu - P22) 
i {Pii + P22) + P12 ~ P21 



(19) 
(20) 



It is instructive to compare the quantum mechanical ap- 
proach and the transfer formalism using an example. Let 
us consider two impurities residing side by side. In this 
instance we know that the transmission and reflection 
amplitude are Eqs. H13|l and H14|) . In the transfer matrix 
formalism, one needs to calculate P = S^^MiMqS to 
obtain the transmission amplitude 



T 



2isin(fc) 



t sin(fc) [{Uo - Ek){Ui ~ Ek) -2] + Uo + Ui- 2Ek - cos(fc)([/o - Ek){Ui - Ek) 

I 



(21) 



Note that T [Eq. lO] is not identical with T [Eq. if^ ]. 
However, while not obvious, they merely differ by a phase 
factor and their magnitudes are the same: T = e^'^^T. 
Clearly the transfer matrix method 'keeps track' of the 
two lattice spacings traversed as the particle is trans- 
mitted to the other side. A similar relation holds for 
the reflection amplitude. Consequently, the transmission 
(reflection) probabilities are identical; |r|^ = |T|^ and 
|i?|'^ = |7?.|'^ in the two formalisms. 

In previous work^^ using the transfer matrix formal- 
ism, we derived a relation between TZ and 7?.', where TZ' 
is the reflection amplitude for the reverse impurity con- 
figuration. It is 



T* {W - TZ) 



2ik 



(22) 



In the example of the two impurities, TZ' is given by TZ 
with Uo and C/i exchanged. Using Eqs. l(T^ and ^^, 
one can show that T and R also satisfy the same relation; 
namely, 



T* {R' - R) 
T{R' -R)* 



2ik 



(23) 



Introducing a phase difference between R and R' such as 
R' — e^R, we obtain 



2ik 



R*T 

RT* 



(24) 



For the case of two impurities we considered, the phase 
difference 5 can be determined by 



tan((5 - 2k) 



2/it^ 
v^ — /i^ 



where 



M = 

V = 



Uo + Ui cos(2fc) + UqUi cos(fc) 
Ui sin(2fc) + UqUi sin(fc) 



(25) 



(26) 
(27) 



Later, we will show that the derivative of S with respect 
to k can be attributed to a phase shift (or time delay) of 
wave functions for the reverse configuration of impurities. 



V. NUMERICAL DIAGONALIZATION AND 
PHASE SHIFT OF WAVE PACKETS 

We have already alluded to the numerical approach for 
a free wave packet transfer in section II. Since the opera- 
tor C^' creates an electron at site i, the initial wave packet 
can be written as |^(0)) = J2i '/'(^^i, 0)Cj^|0), where 



(p{xi,0) 



(27ra2)i/4 



piko{xi-xo) -l{xi-xo)^ /a'^ 



(28) 



As mentioned in Sec. II, xo is the mean position, fco is 
the mean momentum, and a is the initial uncertainty as- 
sociated with the position. If a is much larger than the 



size of the potential region, say /, the wave packet acts 
hke a plane wave when it is scattered off the potential. 
To see the real time propagation of a wave packet in a 
lattice with Nq sites in total, one needs to diagonalize the 
Hamiltonian Eq. lO, which is an (A^o x ^o) raatrix. Since 
the impurity potentials are real, the Hamiltonian matrix 
is real and symmetric. The numerical diagonalization is 
done using the expert driver DSYEVX contained in the 
LAPACK package, which provides either selected eigen- 
values and eigenvectors or the entire spectrum. Using 
the eigenstates \n) and eigenvalues e„ obtained from the 
diagonalization, one can express the wave packet at time 
t as follows: 



Fig.4 (Kim et al) 



No 



|*(t)) = ^|n)(n|vl/(0))e- 



(29) 



The wave packet initially at xq moves to the poten- 
tials and scatters off impurities. In general, the wave 
packet is partially reflected and partially transmitted. 
The mathematical definitions of the reflection and trans- 
mission probabilities are |i?p = J2i<o\'^(''^i'^)\'^ ^^'^ 
|rp = J2i>i Ivi^iit)]"^: respectively as i ^- oo. 

Let us consider the two impurity case again to illus- 
trate the time evolution of a wave packet in the presence 
of impurities in a lattice. The impurity potentials are set 
to be /7o = 1 and C/i = 3 in units of the hopping constant 
to- We consider wave packets with average momentum 
that varies from ko = O.Stt to 0.9tt. The time elapsed 
for the scattering processes to 'finish' depends on ko; for 
example, for fco = O.Gtt it is 160 in our dimensionless time 
unit, while for feg = O.Stt, it would be 240. We consider 
two impurity configurations, [I] and [II]. For [I], we have 
{Uq, Ui) and for [II], {Ui, Uq). As we discussed earlier, the 
reflection amplitude will differ correspondingly; namely, 
R for [I] while R' for [II]. Since T = T' , there is no phase 
shift for the transmitted wave packet as shown in Ref. 
I21I On the other hand, the phase difference S induces 
the phase shift for the reflected wave packets. This can 
be explained as follows: Consider a wave packet ^{x, 0) 
moving with fco to an impurity region 



V'(x,0)= / dfc 5(fc) e^'^o^^-^") 



(30) 



where g{k) ~ ^-"(''-''a) _ After the wave packet scatters 
completely off the impurities at time tg, the wave function 
at t after tg would be 

'il;[i]{x,t) = /"dfc.9K(A:)e-*'=«(^-^«)e~*-^'=* 



dk 5T(A:)e''="(^-^^)e-*^^-* (31) 



where an elastic scattering is assumed, gu ^ R{k)g{k) 
while gT ~ T{k)g(k), and xr{xt) is the mean position 
of the reflected (transmitted) wave packet at tg ■ For the 
reverse configuration [II], 




100 200 

lattice site 



FIG. 4: Phase shift of the reflected wave packets in a lattice 
with two impurities embedded side by side with Uo = 1 and 
C/i = 3 in units of the hopping constant. Initially, the wave 
packet is at xq ~ 150. The two impurities are at 300 and 301. 
The uncertainty parameter a is set to be 20. For |'(/'[/],fl|^ 
(solid curves), the impurity configuration is {Uq, Ui) while for 
!''/'[//], -Rp (dashed curves), it is {Ui,Uo). 



+ dkgT'{k)e'''"^''"''^\ 



lEut 



(32) 



Since R' = R e , the reflected part of ipiu] (a;, t) becomes 

i/'[//],fl(a;,t) -- I dk R{k)e-°'''^~^"^\-'^'>''''~'-'"^-''^''^+'^ 

(33) 
In order to see the phase shift of the reflected wave pack- 
ets, one needs to calculate |V'[/7]._rI^- Expanding E{k) 
and 5(k) around fco and some algebra yields 



\ii]M 



exp 



{x-xr + vpt - dkSp) 



(34) 



where Vq = dkEp. Rigorously speaking, this approxima- 
tion is valid when R{k) ~ R{kQ). Note that the term with 
Eq and i^o disappears when [■(/'[/7]._r|^ is calculated and the 
position of the scattered wave is xr — v^t. Eq. (|34|l in- 
dicates that the phase shift of the reflected wave packets 
is determined by i9fei5o, and the spreading depends not 
only on i9^i?o but also on 9^5o- To be speciflc we define 
a phase shift as the difference between the two reflected 
wave packets for [I] and [II] at their half width as in Ref. 
I2I 

In Fig. 4, we plot the reflected wave packets for [I] 
(solid curves) and [II] (dashed curves) in the case of two 
impurities with Uq = 1 and C/i = 3. The average mo- 
menta for the three wave packets shown are fco = O.Gtt, 
O.Ttt, and O.Stt. As one can see clearly, the phase shift 
between |V'[/],i{|^ (solid curves) and |V'[//],i?P (dashed 
curves) depends on fco. The phase shift for fco = O.Gtt 
is not signiflcant while it becomes bigger as fco increases 



Fig. 5 (Kim et al) 



■u -2 
a 

a 

« -3 



' 1 ' 1 


i 1 i 1 i 


^^* 


X : 


♦ : phase shift 


\ ; 


1 . 1 






0.4 0.6 

k/jt 



FIG. 5: Phase difference (blue dashed curve), its deriva- 
tive (red sohd curve), and the phase shift obtained from the 
wave packet simulations (green diamond symbols) . The phase 
shift is measure by the difference between the reflected wave 
packets t/'[/],Hl s-nd V'li/l.fll ^-t their half width. The impurity 
potentials are Uo ~ 1 and f/i — 3. The phase shift and dS/dk 
agree with one another very well. 



to O.Stt. We also plot the phase difference, its deriva- 
tive, and the phase shift of the reflected wave packets in 
Fig. 5. One may use Eq. (|25l) to obtain the phase dif- 
ference and its derivative in this simple case. However, 
when we consider many impurities we need to use the 
transfer matrix formalism. As shown in the plot, the 
derivative of S and the phase shift obtained numerically 
are in excellent agreement. 



VI. SPIN FLIP SCATTERING 

So far we have discussed scalar potential scattering on 
a lattice. In this section we describe how to study spin flip 
scattering on a lattice. Let us consider the Hamiltonian 

H^^t„ Y. C+Q,-2Jo^arS(-2Ji^SrS,+i 

<i,j>(T I I 

(35) 
where C^ creates an electron with a spin a at a. site 
i, S; is a spin operator located at a site Z, to is a hop- 
ping amplitude between the nearest neighbor sites, and 
Jo is the coupling between an electron and a local spin. 
The electron-spin coupling is assumed here to be a purely 
local (i.e. on-site) interaction. Ji is a (Heisenberg ex- 
change) coupling between two neighbouring spins. This 
model can be used to understand the spin transfer dy- 
namics between an itinerant electron and a ferromagnetic 
spin chain with Ng local spins (5 — 1/2) arranged in a 
one dimensional lattice. To study the spin chain, one can 
extend the quantum mechanical approach as in Ref. Il4l 
or one can extend the transfer matrix formalism^ Here 
we follow the diagonalization method to see the time evo- 



lution of the spins. 

Suppose we send a wave packet representing an elec- 
tron with spin aligned in the +Z direction towards a fer- 
romagnetic spin chain where all spins are aligned in the 
—Z direction. Such a state of the chain will be denoted 
by \G). The incoming wave packet can be constructed as 
follows: 



lV'(0))-5]^(x.,0)C+|0) 



(36) 



Then, the initial wave function of the total system in- 
cluding the incoming electron and the chain is |^(0)) = 
|V'(0))|G). We introduce the total spin operator J = 
a + ^j S/. The Z component, Jz, of the total spin is 
conserved. Hereafter we assume spin 1/2 for both the 
electron and the spins, for simplicity. Since the initial 
value of Jz is (1 — Ns)/2, the possible spin bases would 
be |+)|G) and \—)Si+\G), where Si+ flips the local spin 
at / in the chain. Alternative spin bases could be used by 
utilizing the total spin and its Z component. Including 
the location of the incoming electron, the bases states we 
use are C+^\0)\G) and C+^\0)Si+\G). Now, the Hamilto- 
nian becomes an No{Ns + 1) x No{Ns + 1) matrix. Note 
that the dimension of the Hamiltonian depends only on 
N and Ng and does not depend on the locations of the 
local spins. Even if we include impurities, the dimension 
of the Hamiltonian remains unchanged. To construct the 
Hamiltonian matrix, we need to calculate each compo- 
nent of the matrix. For example. 



(G|(0|Qt {^2Jom ■ SO G+|0)|G) = ^S,,A,i 



Jo 
2 



(G|(0|G,T (-2JiS, . S;+i)G+|0)|G) = -^S.,, 



2 



We therefore solve the eigenvalue problem: H\rjj) — 
Ej\rij), where j = 1, 2, • • ■ No{Ns + 1). Then, we repre- 
sent the time dependent total state using eigenvalues and 



Eit 



eigenstates as follows: \il^{t)) = J2j \Vj){Vj\''P{^))^ 
Alternatively, the total state at t can be expressed as 

No Ns No 

i^w) = E^^w^nio)i^)+EE^MWc,:lio)5,+iG) . 

i=l 1=1 i=l 

.(37) 
Using this expression we can investigate the dynamics of 
a particular spin or the sum of all spins in the chain. For 
example, 

(^(i)|5,,|^(t)) = -iEl^''l'+jEl^''.^'l'(2^M'-l) , 



I'A' 



(38) 
and the total local spin is J^ii^iW)- Since at t = 0, 
V';,i = and J2i Hi? = I7 {Siz{0)) = -1/2 is assumed. 
Thus, (5,(0)) = EiiSiziO)) = -Ng/2 and (J,(0)) = 
{l-Ns)/2. 

By way of an example, consider a chain of three spins, 
ferromagnetically coupled with strength Ji. As an elec- 
tron impinges on the three spin system, they interact 
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FIG. 6: The spin expectation value (Sz) for three spins Si 
{i = 1, 2, 3) and the total spin Stot- The coupling parameters 
are Jo = 2 and Ji — 0. 
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FIG. 7: The spin expectation value (Sz) for three spins Si 
{i = 1, 2, 3) and the total spin Stot- The coupling parameters 
are Jo = 2 and Ji = 1. 



to the same degree, as shown in Fig. 7. Interestingly, 
even though there is an obvious asymmetry (the electron 
strikes the first spin first) the time evolution of the first 
(^i) and the third (S3) spins are almost identical. 
VII. CONCLUSIONS 

In this paper we have described three rather general 
methods to solve scattering problems on a lattice. The 
first two are based on single electron-impurity scattering 
problems, and tell us about the end result of the scatter- 
ing process. We examined several examples of impurity 
scattering. In particular, following Giri et al.S^, we ex- 
plored parameter regions where near resonant transmis- 
sion is expected to occur for the random trimer model. 
As is evident from Fig. 3, there are regions in parameter 
space where near resonant transmission will remain as 
imperfections arise in the trimer itself. It would be very 
interesting to examine the width of the resonance for N 
such trimers arranged randomly on an infinite lattice, as 
N increases. 

The third method developed in this paper is numeri- 
cally intensive, but allows for solution of completely gen- 
eral problems, including those involving spin-flip scatter- 
ing, as explored briefly in the last section. Moreover, one 
can monitor the time-dependence of the wave function 
amplitudes throughout the scattering event. This will 
be of interest as experimental methods evolve to allow 
more controlled temporally and spatially resolved mea- 
surements. We presented, by way of example the case of 
a Heisenberg chain of spins, interacting with an itinerant 
electron (spin current). A natural extension of this cal- 
culation could be to include more spins, so as to model 
an actual magnetized thin film. Generalizing to the case 
of a spin current is straightforward, so long as we assume 
the electrons in the current do not interact with one an- 
other. This will allow us to address detailed questions 
about how the spins in the chain reverse their magneti- 
zation when subjected to a spin current. These and other 
questions will be the subject of future investigation. 



with the electronic spin and change their states. In 
Fig. 6, we first plot the expectation value of Sz for S.; 
(i = 1,2,3), and Stot = J2i=i ^i ^o^ ^^^ uncoupled spin 
case, i.e. with Ji = 0. The electron is coupled to each 
spin with strength Jq = 2, in units of to- Since there is 
no coupling between two nearest spins in this instance, 
each spin evolves independently as a function of time. 
The spin transfer from the incoming electron to the local 
spins occurs mostly for Si while it is minimal for S3. 

On the other hand, when the spins are Heisenberg cou- 
pled with strength Ji = 1 all three spins participate in 
the spin-flip scattering with the incoming electron almost 
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